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Abstract. We present a single-freeze-out model with thermal and geometric 
parameters dependent on the position within the fireball and use it to describe 
the rapidity and transverse-momentum spectra of pions, kaons, protons, and 
antiprotons measured at RHIC at v^sjvjv = 200 GeV by BRAHMS. THERMINATOR 
is used to perform the necessary simulation, which includes all resonance decays. 
The result of the fit to the data is the expected growth of the baryon and strange 
chemical potentials with the spatial rapidity oy . The value of the baryon chemical 
potential at q:|| ~ 3 is about 200 MeV, i.e. lies in the range of the highest SPS 
energies. The chosen geometry of the fireball has a decreasing transverse size as 
the magnitude of oy is increased, which also corresponds to decreasing transverse 
flow. The strange chemical potential obtained from the fit to the /K~ ratio 
is such that the local strangeness density in the fireball is compatible with zero. 
The resulting rapidity spectra of net protons are described qualitatively within the 
statistical approach. As a result of our study, the knowledge of the "topography" 
of the fireball is acquired, allowing for other analyses and predictions. 
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More details of the presented material can be found in Ref. [1] . 

An important goal of relativistic heavy-ion physics is to find the "topography" of 
the system at freeze-out, i.e. to determine the dependence of thermal parameters on 
the location, in particular on the spatial rapidity Q!|| = ^ log {t — z)/{t + z). Numerous 
studies of the global particle ratios [2-18] can be divided into two categories: the An 
studies at low energies (SIS, AGS) and the analyses at mid-rapidity for approximately 
boost-invariant systems at high energies (RHIC). In the latter case the yields dN/dy 
vary only by a few percent in the window \y\ < 1. When the system is not boost- 
invariant the task is more involved because particles detected at a given pseudorapidity 
may originate from different locations in the fireball. Thermal conditions and flow 
change from place to place, which should be taken into account. In addition, the 
resonance decays are more complicated to treat. As a result, a full-fledged simulation 
is necessary. We use THERMINATOR [19] with a suitably extended single-freeze-out 
model of Ref. [11,12,20]. We recall the approach reproduces in an efficient way 
the transverse- momentum spectra, including strange particles [20], the production 
of resonances [21], the charge balance functions [22], the elliptic flow [23], the HBT 
radii [24], and the transverse energy [25,26]. To fix the model parameters we use the 
data from the BRAHMS collaboration at y^s^N = 200 GeV [27-30]. 

Our extension to boost-non-invariant systems consists of two elements: the choice 
of the freeze-out hypersurface and the collective flow, as well as incorporation of the 
dependence of thermal parameters on the location. We use the freeze-out hypersurface 
parameterized as 

— T (coshaxcoshaii , sinhaj^ cos (/>, sinhax sin 0, cosha^sinhaii) . (1) 

The parameter r is the proper time at freeze-out, and is related to the transverse 
radius, p = \J ^y"^ = rsinhaj^. The four- velocity is chosen to follow the Hubble 
flow = /t. In the boost-invariant model p was limited by the space- independent 
parameter pmax, or < oi±_ < a™^'^. Now we take 

< < ar'^(aii) ^ a--(0)exp (^^^^ • (2) 

The interpretation of Eq. ([2|) is following: as we depart from the center, increasing 
)q!|| ], we simultaneously reduce a±, or pmax- The rate of this reduction is controlled by 
a new model parameter, A. In other words, the fireball gets thinner and thinner as we 
depart from the central rapidity region. The geometry and expansion of the fireball are 
described by three parameters: r, controlling the overall abundance of particles, pSax, 
influencing the slope of the pT-spectra, and A, controlling the fall-off of the rapidity 
spectra. With the standard parameterization of the particle four- momentum in terms 
of rapidity and the transverse mass, p'* = (r7i^cosh?;,p^ cos ip,p± simp, mj^sinhy) , we 
find 

p ■ u — m^cosh(Q;^)cosh(a|| ~ y) — p_Lsinh(a_L) cos(0 — ip), 

d'^I] ■ p ~ T^da^id(f>smha±da±p ■ u (3) 

where d^'S,^ is the volume element of the hypersurface. These expressions are used the 
Cooper-Frye [31] formula for the momentum density of a given species of primordial 
particles: 

roo ^q^'="'(q||) ^27r 



2npTdpT dy 



day / da±_ d(j)p ■ uf {Pp ■ u - (3p{a\\)) , (4) 
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Figure 1. Left; the model baryon and strange chemical potentials plotted as 
functions of «|| . Parameters of Eq. JSJ are obtained from the fit to the BRAHMS 
data [27,28]. The points represent a naive calculation based of Eq. (|6)l. Right: 
Comparison of obtained from the fit to the data (solid line) and from /ig and 
the condition of zero local strangeness density (dashed line). 



where p ■ u is taken at = (azimuthal symmetry), / denotes the 
Bose- Einstein or Fermi-Dirac statistical distribution function, /3 = 1/T, and 
/Lt(a||) = B/iB(Q!||) + S'/xs(a||) + /3/j,/3(q!||), with S, and I3 denoting the baryon 
number, strangeness, and the third component of isospin of the given particle. Thus 
we admit the dependence of chemical potentials on ay. In general, the temperature T 
also depends on a|| . In this work, however, we apply the model for not too large values 
of rapidity, \y\ < 3.3, where the obtained values for /xb are less than ^ 250 MeV. The 
universal freeze-out curve [8] gives a practically constant value of T in this range, thus 
we fix r = 165 MeV everywhere. 

The functional form of the chemical potentials is parameterized as follows: 

=Ai,(0) + , i^B,S,h. (5) 

Practice shows that to a good approximation the statistical distributions may be 
replaced by the Boltzmann factors. Then the integrand of Eq. ^ contains the term 
exp[— /?mj^cosh(Q;j_)sinh(Q;|| — y) + /3/i(Q;||)]. The relevant integration range in a|| is 
sharply peaked around ay ~ y, thus the chemical potentials are taken approximately 
at /ii(a||) ~ and the factors exp[/3/x(?/)] can be pulled out from the integration. 

As a result, the following approximate relations for the ratios of yields hold: 

p K+ 
- ~ exp(2/3^i3), — 
p K 

Upon inversion usiy) — ■^T\og{p/p), etc. With the help of these approximate 
equalities we set the starting values of the parameters /ii(0) and Ai, which are then 
iterated by THERMINATOR to fit the data. The A parameter is fixed with the pion 
rapidity spectra dN^±/dy, with the optimum value of A = 3.33. The other geometric 
parameters are taken from earlier fits to mid-rapidity p^-spectra: r = 9.74 fm, 
Pmax(z = 0) = 7.74 fm. 

The result of our optimization is shown in Fig. [1] The optimum parameters are: 

fisiO) 19 MeV, ^is{0) = 4.8 MeV, fii^iO) = -1 MeV, 
Ab = 0.65, As = 0.70, Ai^ = 0. (7) 



exp(2/3^i3), ~ exp (2/3^5), — ~ exp( 2/3^/3). (6) 
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Figure 2. Dependence of particle ratios on rapidity: left - p/p, right — /K~ . 
The open triangles are the Brahms data [27,28], while the filled circles show the 
result of the model simulation. The data for p and p are not corrected for the 
feed-down from weak-decays [27]. 



We observe the expected behavior for the baryon chemical potential, which increases 
with |a|||. The value at the origin is 19 MeV, somewhat lower than the earlier mid- 
rapidity fits made in boost-invariant models in Refs. [11,32], yielding 26 JVIeV. The 
lower value in our case is well understood. The previous mid-rapidity fits include the 
data in the range \y\ < 1. This range collects the particles emitted from the fireball at 
I a II I < 2, hence the value oi fj,B in the previous mid-rapidity fits is an average of our 
fiB{c(\\) over the range, approximately, |q;||| < 2, with some weight proportional to the 
particle abundance. A similar effect occurs for /zg. We do not incorporate corrections 
for the feed-down from weak decays (i.e. all the decays are included), since this is the 
policy of Ref. [27] for the treatment of p and p. 

We note that at a|| = 3 the value of /is is 200 MeV. This value is comparable to 
the highest-energy SPS fit (y/s^N = 17 GeV), where fiB ^ 230 MeV. The behavior 
of the strange chemical potential is qualitatively similar. It also increases with |a|||, 
growing from 5 MeV at the origin to 50 MeV at a\\ = 3. The ratio piB{a\\)/ iJ,s{o!.\\) is 
very close to a constant ~ 4 — 3.5. The dots in the left panel of Fig. [T] show the result 
of the naive calculation of Eq. ([6]). We note that these points are close to the result 
of the full-fledged fit of our model. This is of practical use, since the application of 
Eq. ([6]) involves no effort (as done e.g. in Ref. [33-35]), while the full model simulation 
incorporating resonance decays, fiow, etc., is costly. 

There is another interesting point. In thermal models one may obtain the local 
value of the strange chemical potential, f^g, at a given fiB with the condition of 
the vanishing strangeness density, ps = 0. The result is shown in Fig. [Tl where we 
compare the strange chemical potential obtained from the fit to the data (solid line) 
and from the condition of zero local strangeness density at a given /^^(aii) (dashed 
line). The two curves turn out to be virtually the same. This shows that the net 
strangeness density in our fireball is, within uncertainties of parameters, compatible 
with 0. This is not obvious from the outset, as the condition of zero strangeness density 
is not assumed in our fitting procedure. Although this feature is natural in particle 
production mechanisms, in principle only the total strangeness, integrated over the 
whole fireball, must be initially zero. Variation of the strangeness density with a|| is 
admissible, but turns out not to occur. 

Figure [5] shows the quality of our fit of the parameters of the chemical potentials, 
Eq.®, for the ratios p/p and / . In Fig. [3] we compare the obtained rapidity 
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Figures. Left: Left: Rapidity spectra of 7r+ , /f"*", and . Tlie data come from 
Refs. [27,28] (circles - 7r+, squares - , triangles — K~), while the histogram 
lines show the model results. For 7r+ the solid (dashed) line corresponds to the 
full feeding (no feeding) from the weak hyperon decays. The experimental pion 
yields are corrected for weak decays [28]. Right: spectrum of net protons, p — p. 
The data points come from Refs. [27, 28] , while the solid (dashed) histogram lines 
show the result of the model simulation with full feeding (no feeding) from the 
weak hyperon decays. Data points should be compared to the model with full 
feeding (solid lines). 

spectra of tt^, -ftT^, and to the experiment. The experimental yields for the pions 
are corrected for the feed-down from the weak decays as described in Ref. [28]. For 
that reason for the case of tt"*" we give the model predictions with the full feeding 
from the weak decays (solid line) and with no feeding at all (dashed line). We note 
a reasonable agreement, with the data falling between the two extreme cases. The 
spectra of kaons are also reproduced reasonably well. We see from the right panel of 
Fig. [3] that the qualitative growing of the net-proton spectrum with y is reproduced. 
Note that predictive power is left for this observable, as only the ratio p/p is used to 
&t hb- The model curves with and without the weak-decays correction embrace the 
experimental points. 

Here are the highlights this talk: 

(i) Decreasing yields of particles with rapidity require that the system becomes 
colder, thinner, or more dilute as one departs from mid-rapidity. In our analysis 
for the highest RHIC energies we take it to be thinner, keeping the temperature 
constant. 

(ii) The naive extraction of the baryon and strange chemical potentials from ratios 
of p/p and / according to Eq. ([6]) works surprisingly well, as shown in the 
comparison to the full calculation in Fig.[T] This is because at ^/snn — 200 GeV 
the dependence of yields on a\\ does not depart much from linearity within one 
unit of rapidity, which is the window contributing to the a\\ integration in Eq. 

In addition, the feeding contributions from resonance decays to both particles in 
the ratio are approximately proportional to each other. At lower energies (SPS) 
the naive extraction should not work. 

(iii) The baryon and strange chemical potentials grow with a||, reaching at y 3 
values close to those of the highest SPS energies of ^/snn = 17 GeV. This agrees 
with the conclusions of Roehrich [36] . 

(iv) At mid-rapidity the values of the chemical potentials are even lower than derived 
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from the previous thermal fits to the data for \y\ < 1, with our values taking 
Ij,b{0) = 19 MeV and /i5(0) = 5 MeV. The reason for this effect is that the 
particle with \y\ < 1 originate from a region |a||| < 2, and on the average the 
effective values of chemical potentials are larger compared to the values at the 
very origin (cf. Fig. [T]) . 

(v) The local strangeness density of the fireball is compatible with zero at all values 
of q;||. Although this feature is natural in particle production mechanisms, here it 
has been obtained independently just from fitting the chemical potentials to the 
data. 

(vi) The ratio of the baryon to strange chemical potentials varies very weakly with 
rapidity, ranging from ~ 4 at midrapidity to ^ 3.5 at larger rapidities. 

(vii) The d'^N/{2Trp±dp±dy) spectra of pions and kaons are also well reproduced [1], 
complying to our hypothesis for the shape of the fireball in the longitudinal 
direction. 

(viii) The feature of the increasing yield of the net protons, p — p, with rapidity is 
obtained naturally, explaining qualitatively the shape of the rapidity dependence 
on purely statistical grounds. 
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